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We present a new approach to quintessential inflation, in which both dark energy and inflation 
are explained by the evolution of a single scalar field. We start from a simple scalar potential 
with both oscillatory and exponential behavior. We employ the conventional reheating mech- 
anism of new inflation, in which the scalar decays to light fermions with a decay width that 
is proportional to the scalar mass. Because our scalar mass is proportional to the Hubble 
rate, this gives adequate reheating at early times while shutting off at late times to preserve 
quintessence and satisfy nucleosynthesis constraints. 



1 Intoduction 

If we assume the correctness of the standard Friedmann equation evolution, the existence of 
dark energy engenders two profound dilemmas. The first is the cosmological constant problem, 
and the second is the "why now?" problem. Quintessence models attempt to address the second 
problem by introducing a very weakly coupled scalar field whose potential and/or kinetic function 
have special properties. One of the most successful approaches to quintessenc e^^ is to combine 
tracking with an oscillating behavior in the quintessence potential. In such models the equation 
of state parameter w(z) has a periodic component, leading to occasional periods of accelerated 
expansion during epochs where w(z) ~ —1. 

It is natural in this context to ask whether the quintessence scalar could replace the infla- 
ton. The idea of quintessential inflation has been examined by a number of authors ^^ l The 
straightforward approach is to cobble together a scalar potential which has both a flat, large 
vev portion (for inflation) and a flat, small vev portion (for quintessence). These features are 
conn ected by a steep step which corresponds to a period of cosmic kination. As discussed in 
I3H3 such models suffer from generic problems. First, they require significant ad hoc tuning 
to simultaneously produce the features of inflation and quintessence. Second, they require a 
"sterile" inflaton, in order to avoid the decay of the putative quintessence scalar at the end of 
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Figure 1: The upper (lower) panel shows the scalar potential V{9) (V(4>)) for 6=1, plotted on a logarithmic 

scale. 



inflation. This in turn requires new mechanisms for reheating, such as gravitational de Sitter 
phase particle production, leading to difficulties in satisfying the constraints of CMB anisotropics 
and of primordial nucleosynthesis. 
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Our approach to quintessential inflation 1 11 is to take advantage of the tracking and oscil- 
latory potential features that work so well in addressing the "why now" problem of quintessence 
alone — . We will describe a model with a simple scalar lagrangian with exponential and oscil- 
latory features. The model uses the conventional reheating mechanism of new inflation!^, in 
which the scalar decays to light fermions with a decay width that is proportional to the scalar 
mass. We show that the scalar mass is proportional to the Hubble rate. As a result, the model 
has adequate reheating at early times while naturally shutting off at later times. 

We present a successful model with only three adjustable parameters. One parameter con- 
trols the period between inflationary epochs, a second controls the overall decay width, and the 
third parametrizes our ignorance about the relative fraction of matter versus radiation produced 
by reheating. These three parameters are adjusted to produce sufficient inflation along with the 
correct fractions Q r /ri\, ^m/^A of radiation, matter, and dark energy, as measured today. 
Having thus fixed the model we find that we automatically satisfy all constraints of primordial 
nucleosynthesis, CMB, and large-scale structure. 



2 An oscillatory potential 



We start with a simple model with a single real scalar quintessence field 9. The action is 



\f{e)sTd^ed v 6 - v(9) 



(i) 



where the kinetic function f(9) and potential V(9) are given by: 
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where Mp is the Planck mass: 1.22 xlO 19 GeV; po is the dark energy density observed today: 
~ (10~ 4 eV) 4 ; b is a dimensionless parameter which controls the periodic behavior. A canonical 
kinetic term can be restored via a field redefinition 9(x) — ► (j)(x), where 



<p(x) = 4>o cos 9 , 



(4) 



with (j) Q = ^2,M 2 /irb 2 . 

In the approximation where we ignore the energy density of radiation and matter, and 
where the only friction is from the metric expansion, the evolution of the model can be solved 
analytically. Energy conservation requires: 
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where H is the Hubble rate, w is the equation of state parameter, and pg is the dark energy 
density. The solution to this equation as a function of the scale factor a(t) is: 



p e (a) = p exp 



-3 f -(1 +w(a)) 
Ji a 



(6) 



where po is the dark energy density at a = 1 (today). 
We also know that 



V(6) = -(l-w)p e . 



Making the Ansatz 
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one immediately gets a solution to the (flat) Friedmann equation combined with the relations 
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The expectation value of the quintessence field 9 evolves logarithmically with scale factor 
from a positive initial value to zero today. Accelerated expansion corresponds to epochs (such as 
today) where 9 is evolving through one of the flat "steps" of the potential. From @ we see that 
the kinetic function is simultaneously suppressed in this epochs, slowing the roll of the scalar 
field evolution. 

The equation of stat e param eter w(a) has t he s ame periodic form assumed in recent phe- 
nomenological analyses ^^^^^^ . The analysis in ED showed that, for 6 = 1, w(a) is consistent 
with all current data from observations of the CMB, Type IA supernovae, and large scale struc- 
ture. It follows a fortiori that our model with any choice of b less than one also agrees with this 
data. 

To complete the model, we will assume that the quintessence field (j) has a weak perturbative 
coupling to light fermions. This is the standard reheating mechanism of new inflation 14 . To 
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Figure 2: History of a slinky universe. Shown are the relative energy density fractions in radiation (green), matter 
(blue), and dark energy (red), as a function of the logarithm of the scale factor a(t). 

avoid the strong constraints on long-range forces mediated by quintessence scalars^l, it is sim- 
plest to imagine that our scalar only has a direct coupling to a sterile ne utrin o. This is sufficient 
to hide the quintessence force from Standard Model nonsinglet particles'^, while still allowing 
the generation of a radiative thermal bath of Standard Model particles from quintessence decay. 

With this assumption for reheating the evolution equations for quintessence, radiation, and 
matter become: 



pe = -3H(1 + w)pe - fcom^l + w)p e ; 
p r = -4:Hp r + (l-f m )kom^(l + w)p e ; 
Pm = -3Hp m + fmkom^l + w)p e ; 



(11) 



where ko and f m are small dimensionless constants. As long as 9 is not near a multiple of vr/2, 
it is a reasonable approximation to make the replacement k^m^ — > kH where k is another 
small dimensionless parameter. This replacement decouples the 9 evolution equation from the 
Friedmann equation, giving an immediate analytic solution: 



~(3 + k) (29 -sin 29) 
b 



p g (a) = p exp 

We have used this approximation in the solutions quoted below 



(12) 



3 Results 

Figure [2] shows the results obtained from our model with b = 1/7, k = 0.06, and f m = 10~ n . 
Shown are the relative energy density fractions in dark energy, radiation, and matter, as a 
function of logo. We have chosen to integrate the evolution equations starting from a = 10 -42 , 
which in our model corresponds to a temperature of slightly less than 10 16 GeV, and an initial 
comoving Hubble radius of about 100 Planck lengths. For simplicity we have also chosen the 
value of w(a) now to be exactly —1. Neither of these choices corresponds to a necessary tuning. 

Our three adjustable parameters have been chosen such that the values of Q r /U\, £l m /VL\ 
come out to their measured values at a = 1, and such that we have sufficient inflation. The 
latter is checked by computing the ratio of the fully inflated size of the intitial comoving Hubble 
radius to the current comoving Hubble radius. This ratio is about 3 in our model, indicating 
that the total amount of inflation is indeed enough to solve the horizon problem. The flatness 
problem is solved because the total Q r + 0, m + £7a = 1 within errors. 

From Figure (J2J) we see that we are currently beginning the third epoch of accelerated 
expansion. The first epoch of inflation accumulated about 18 e-foldings. Quantum fluctuations 
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Table 1: Relative density fractions of dark energy, radiation, and matter, as a function of the scale factor. Also 
shown are the temperature T in GeV, and the equation of state parameter w(a). 

during this epoch produced the spatial inhomogeneities responsible for large scale structure and 
CMB anisotropics observed today. Constraints on the physics responsible for the primodial 
power spectrum of these density fluctuations can be set with WMAP and 2dF data, under the 
assumption that the Hubble rate is dominated by the contribution from oa, during the observable 
part of inflation^. As can be seen from the fignre a second period of accelerated expansion 
began just before the electroweak phase transition, and ended well before BBN. The temperature 
history near the EWPT is shown in Figure 

Also shown is the Hubble parameter H of the model normalized to the expansion rate H ra< i 
for pure radiation. H ra d corresponds to what is assumed in the standard paradigm. Notice that 
for temperatures of a few GeV the expansion rate is actually somewhat larger than normal, but 
at higher temperatures it is much less than normal. 

Such a nonstandard thermal history will impact on electroweak baryogenesis. For a Higgs 
sector such that the EWPT is first order, the change in the net baryon asymmetry is proportional 
to -\og(H / H ra d) , where H is the expansion rat e d uring the phase transition, and H ra d is the 
corresponding expansion rate for pure radiation 1 ^. If the Higgs sector is such that the EWPT 
is second order, the baryon asymmetry is proportional to the expansion rate^l Clearly one 
should revaluate the popular scenarios for electroweak baryogenesis in this light. 

Such a model will have major implications for predictions of the relic abundance of dark 
matter particles with Terascale masses. The dominant production mechanism for such particles 
may be scalar decays, as suggested by Figure Even if the dark matter particles are thermal 
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Figure 3: The temperature history of Model 2 (red/solid) near the electroweak phase transition. Shown in 
green/dashed is the Hubble parameter H of Model 2 normalized to the expansion rate H ra d for pure radiation. 

relics, their abundance now will be affected by the nonstandard expansion rates at earlier times. 

To compare with contraints from primordial nucleosynthesis-™, we note that the second 
epoch of dark energy domination ended well before a ~ 10 -10 , the time at which nucleosynthesis 
occured. Indeed dark energy reheating effects are completely negligible from a ~ 10~ 10 until 
today. 

Figure © also shows one (very brief) prior epoch of matter domination. This only occurs 
if there is some suitably heavy and long-lived matter around to go out of thermal equilibrium 
when a ~ 10 -22 , e.g. a superheavy neutrino. 

It is difficult to extract a precise prediction for the spectral indices of this model, since we 
are never strictly in the slow roll regime. We will be content here with a rough estimate. This 
is obtained starting from a canonical field redefinition: 



where H{(ft) is the Hubble rate we would get ignoring radiation and matter. During inflation, H 
is approximately constant, but this is not an especially good approximation since we are not in 
a slow roll regime. This is similar to the oscillatory models of quintessential inflation discussed 
in I2SI. Taking the potential in the inflationary phase to be approximated by V tx (ft 2 , we can 
estimate the scalar spectral index n s : 




(13) 



In terms of the canonical scalar <ft, the potential can be written: 



V{<ft) oc (ft 2 H 2 (<ft), 



(14) 



(15) 



with iV the number of e-folds between the Hubble radius exit and the end of the second in- 
flationary period. For our model N = 29, giving n s = .93, in good agreement with recent 
observations. 

It is fair to say that, compared to the simple model presented above, the standard new 
inflation scenario looks rather extreme. In the evolution history portrayed in Figure ^ the 
interplay between inflationary expansion and reheating is much milder. In fact, apart from a 
few very brief periods, reheating effects in our model do not actually increase the temperature 
of the thermal radiation bath. Instead, the temperature is almost always decreasing, but it 
decreases more slowly than in the standard ACDM evolution. 

Our particular realization of Slinky inflation must be regarded with some care, since <f>o > Mp 
(as in all models of chaotic inflation) and we have not invoked any consistent Planck scale frame- 
work such as string theory. However the model itself is surprisingly simple, and the physical 
picture which emerges from it has some compelling features. These are worthy of further inves- 
tigation. Also, since the target space parametrized by our scalar is S 1 , it would be interesting 
to extend this scenario to a class of nonlinear sigma models with other compact target spaces. 
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